Spin excitations and mechanisms of superconductivity in cuprates 
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A microscopic theory of spin excitations in strongly-correlated electronic systems within the t-J 
model is discussed. An exact representation for the dynamic spin susceptibility is derived. In the 
normal state, the excitation spectrum reveals a crossover from spin-wave-like excitations at low 
doping to overdamped paramagnons above the optimal doping. At low temperatures, the resonance 
mode at the antiferromagnetic wave vector Q = 7r(l, 1) emerges which is explained by a strong 
suppression of the spin excitation damping caused by a spin gap at Q rather than by opening of 
a superconducting gap. A major role of spin excitations in the d-wave superconducting pairing 
in cuprates is stressed in discussing mechanisms of high- Tc superconductivity within the Hubbard 
model in the limit of strong correlations, while electron-phonon interaction and a well-screened weak 
Coulomb interaction are not essential. 
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I. INTRODUCTION 

Recent studies of electron and spin-excitation spectra 
using angle-resolved photoemission (ARPES) and inelas- 
tic neutron scattering (INS) have revealed an important 
role of antiferromagnetic (AF) spin excitations in the 
|"kink" phenomenon and the d-wave pairing in cuprates. 
In particular, in Ref. [l| quantitative analysis of the 
AF spin-excitation spectrum measured by INS and of 
ARPES data for the spin-fermion coupling of the same 
YBa2Cu306.6 (YBCOg.e) crystal were used for numeri- 
cal solution of the Eliashberg-type equations. The super- 
conducting transition temperature was found to exceed 
% = 150 K. 

' The main argument against the spin-fluctuation pair- 
ing mechanism, a weak intensity of spin fluctuations at 
the optimal doping seen in INS experiments 0] , was dis- 
missed in recent resonant inelastic x-ray scattering Q. 
In a large family of cuprates paramagnon AF excitations 
with dispersions and spectral weights similar to those of 
magnons in undoped cuprates were found. A numer- 
ical solution of the Eliashberg equations for the mag- 
netic spectrum found in YBCO7 and for the electron- 
spin interaction described by the t-J model results in 
% = 100 - 200 K. These calculations based on experi- 
mental data demonstrate that spin fluctuations have suf- 
ficient strength to mediate high-temperature supercon- 
ductivity in cuprates and, therefore, alternative mech- 
anism based on electron-phonon interaction (EPI) (see, 
e.g., Q) seems to play a secondary role in cuprate super- 
conductivity. Strong EPI observed in polaronic effects in 
cuprates may be irrelevant for the d-wave pairing medi- 
ated by / = 2 component of EPI as pointed out in Ref. ^ . 

In this report we briefly consider a microscopic theory 
of spin-excitation spectrum in strongly correlated elec- 
tronic systems (SCES) [1, 0] ■ Using a model for the spin- 



excitation spectrum, we consider spin-fluctuation pairing 
within the Hubbard model in the limit of strong corre- 
lations, > t 0, i]. To compare various mechanisms 
of superconducting d-wave pairing, we take into account 
also EPI and a well-screened weak Coulomb interaction 
considered in Ref. [l3l. We show that the latter gives 
a small contribution for the d-wave pairing and cannot 
suppress the superconductivity. 



II. SPIN-EXCITATION SPECTRUM 

To describe the low-energy spin excitations in SCES 
the one-subband t-J model can be used: 
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where is the hopping integral and is the exchange 



interaction. Here cL 



(1 



-a) are the pro- 



jected Fermi operators acting in the the singly occu- 



Sf 



pied subband and Ui = J2a'"'i,<^^ ''^i-' 
(1/2) cJ^T"^,CiCT' are the spin-1/2 operators where 
t", are the Pauli matrices, a — ±1. 

Using the projection technique for the Kubo-Mori re- 
laxation functions, an exact representation for the dy- 
namical spin susceptibility (DSS), the retarded Green 
function (GF) of the transverse spin-density operators 

(see also (llj): 



S'q = S'q ± iS^ , can be derived 



x(q,c.) = -((5+|5lq)). 



m(q) 



(2) 



where m(q) = {[iS+,SZ^]} - {[[S+ , H], SZ^]). The 
static spin-excitation spectrum ujq is calculated from 
the equality for Kubo-Mori correlation function m(q) — 



. ^ . >, (S'+,S'_q), where (-5*+, S-.q) is evalu- 

ated in a generalized mean- field approximation t6!|. The 
self-energy is given by the retarded GF, 
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FIG. 1: Spectral function at 5 = 0.2 compared with experi- 
mental data ^] at T = 5K (squares) and T — lOOK (circles). 
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FIG. 2: Spectral function at 5 = 0.09 compared with experi- 
mental data 12\ at T = 8K (squares) and T — 85K (circles). 



The "proper part" (pp) of the GF(I3]) describes the pro- 
jected time evolution as in the original Mori projection 
technique. The self-energy ([3]) is defined in terms of the 
force operators -Sf = [[Sf, {Ht + Hj)], {Ht + Hj)] = 
X^a ^i' ("^ — ■^"^)' where Ht and Hj are the 

hopping and the exchange parts of the Hamiltonian ([1]). 

In the Heisenberg limit at zero doping, (5 = 0, the 
self-energy is determined by the force F-'^ . At a fi- 
nite hole doping, 6 > 0.05, the largest contribution 
to the self-energy ^ is given by the hopping term 

Pr = j:,^nU,{t,n [Hr^n+H+,] ~ {^ ^ j)} , where 

H^jn = c\„SjCna + 4^(1 ^ nj^^„)cn-\ ■ Wc Calculate the 
self-energy in the mode-coupling approximation (MCA), 
{clS-Cr..\cl,{t)S+m',{t)) = {clc,.{t)) {S-S+{t)) 
{cnacl^,^{t)) . In the superconducting state, the anoma- 
lous correlation functions {cl _^cl^,^{t)) {S^Sj^{t)) 
{cna-Ci' .-a-{t)) are also taken into account. Using the 
spectral representation for these two-time correlation 
functions both the real, E'(q, w), and the imaginary, 
I]"(q,a;), parts of the self-energy ([3]) are calculated Q- 

The spectrum of spin excitations ujq and the damp- 
ing Fq = — (l/2)I]"(q, Wq) are calculated in a broad 
region of temperature and doping. In the Heisenberg 
limit at ^ = the spectrum of spin excitations reveals 
well-defined quasiparticles with Fq ^ ujq characteristic 
to the Heisenberg model. However, for non-zero dop- 
ing the spin-electron scattering contribution I]"(q, w) in- 
creases rapidly with doping and temperature and already 
at moderate hole concentration far exceeds the spin-spin 
scattering contribution Ej(q, w) . We conclude, that at 
low enough doping and low temperatures well-defined 
spin-wave-like excitations propagating on the AF short- 
range order background are observed, while for higher 
doping and temperatures a crossover to AF paramagnon- 
like spin excitations occurs as found in INS experiments. 

In the superconducting state the spectral function 
x"(Q,w) = Imx(Q,w) were calculated assuming the 
d-wave gap function Aq = (A/2)(cosg2; — cosg^) 0. 
The DSS reveals a pronounced resonance mode (RM) 
at low temperatures due to a strong suppression of the 



damping of spin excitations. This is explained by an 
involvement of a spin excitation in the decay process de- 
scribed by creation of three excitations: particle-hole pair 
with energies uJi + uj2 and a spin excitation with energy 
013 which is controlled by the energy and momentum con- 
servation laws, uj = LUi + UJ2 + and q = qi + q2 + qs- 
Due to the spin gap in the spin-excitation spectrum at 
Q the spin excitation with the energy ^3 ~ wq in this 
process plays a dominant role in limiting the decay of 
the RM in comparison with the superconducting gap in 
the particle-hole excitation. Since luq shows a weak tem- 
perature dependence at T < Tc the RM does not re- 
veal an appreciable temperature dependence and can be 
observed even above Tc in the underdoped region (see, 
e.g., ^M)- 

Figure [T] shows the temperature dependence of the 
spectral functions in the overdoped case at S = 0.2 and 
experimental data (symbols) for YBCO6.92 0- The RM 
having a high intensity at low temperatures strongly de- 
creases with temperature and becomes very broad at 
T ^ Tc- The spectral function for the underdoped 
case S = 0.09 is plotted Fig. H The RM shows a 
weak temperature dependence and is still visible even 
at T = 85 K= 1.4 Tc as found in YBCOg.s crystal fB\. 

Thus, as compared with the spin-exciton scenario for 
the RM based on the random-phase approximation where 
only electron-hole bubble diagrams are taken into ac- 
count (see, e.g., [HI), we propose an alternative expla- 
nation of the RM which is driven by the spin gap at Q 
rather than by opening of the superconducting gap. 



III. SPIN-FLUCTUATION d-WAVE PAIRING 

Despite of intensive search for the mechanism of high- 
temperature superconductivity in cuprates, there is still 
no commonly accepted theory (for a review see 14]). 
A microscopic theory of superconducting d-wave pair- 
ing mediated by AF exchange interaction and spin- 
fiuctuations induced by kinematic interaction has been 
developed within the t-J model in Ref. [l^ and the Hub- 
bard model in Ref. [8j. 
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Recently, the problem of superconductivity in the re- 
pulsive Hubbard model in the weak correlation limit was 
discussed. In Ref. an asymptotically exact solution 
for the d-wave pairing was found, while consideration of 
a well-screened weak Coulomb interaction (CI) has not 
shown a possibility for superconducting pairing [lol |. To 
resolve this controversy, we have considered superconduc- 
tivity in the Hubbard model in the limit of strong corre- 
lations, U ^ t, taking into account also a well-screened 
weak CI and EPI: 
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We introduced here the Hubbard operators (HOs) X"^ = 
|ia)(i/3| for the four states on the lattice site i: an empty 
state a — |0), a one-hole state a = \a) with the spin 
a = ±1/2 — (t,i), — —a, and a two-hole state |2) = 
t4-)- To apply the model for cuprate superconductors, 
we introduce the single-particle energy £i = Sd — ^J- a.s 
an energy of the one-hole d-state. The two-hole energy 
£2 = 2ei + [/ is an energy of the p-d singlet state where 



U 



€d is the charge-transfer energy between the 



oxygen p and copper d states. 

The last term in Q denotes a weak screened CI V{ij) 
between charge carriers in the plane and EPI g{ij) for 
charge carriers 



(5) 



where Uj is a displacement for a particular phonon mode. 
Ni = X^'^ +2XP is the number operator. The chem- 
ical potential /i depends on the average hole occupation 
number n — 1 + 6 = { Ni). 

Using the projection technique in the equation of mo- 
tion method for the GF in terms of the HOs as described 
in Q we can derive an exact Dyson equations for the 
two-subband matrix GFs. The normal GF can be written 



as, 



GA,(k,c^) = c^fo-£(k)-E(k,c^; 



(6) 



where e(k) is the hole energy in the mean-field approxi- 
mation (MFA) and S(k, w) is the normal self-energy. The 
anomalous (pair) GF reads, 

X^^ 



F^iKc) = {{^2^, j \{XZiX%))) 
= -GAr(k, -w) (^cr(k,a;) Gcr(k, w). 



(7) 



The superconducting gap function 1^0- (k,w) — ^aO<-) + 
$o-(k, w) has a nonretarded contribution Acr(k) deter- 



mined by the AF exchange interaction and CI in MFA 
and the anomalous self-energy $cr(k, w). 

The self-energies E(k, w), $cr(k, w) are calculated in 
the MCA by assuming an independent propagation of 
Fermi-like and Bose-like excitations in multiparticle GFs. 
Below we consider the hole-doped case, n = 1 + S > 1. 
The diagonal components of the self-energies for the two- 
hole subband can be written as 

+00 

E22(k,c.) - 1^ 1 dzi^(+)(c.,z|q,k-q) 



[-(lAQ2)ImG22(q,z)], 



(8) 



+ 00 



x[-il/7rQ2)lmF^'{q,z) 



(9) 



where Q2 = n/2 is the weight of the second subband. 
The kernel of these integral equations has a form, similar 
to the strong-coupling Eliashberg theory \v^ : 



+00 



dio' tanh + coth ^ 



(jj ~ z — UJ' 



i^W(c.,.lq,k-q)=| ^ 

||i(q)|2lmx./(k - q,uj') ± l^k-ql^Im Xp/i(k - q,uj') 
± [l^k-qP + |t(q)lV4] Imxc/(k - q,^')}, (10) 



where the spectral density of bosonic excitations are de 
termined by the dynamic susceptibility for spin fluctu- 
ations, Xsfi^,^) — — ((Sq|S_q))(j, charge fluctuations 
Xc/(q,w) = -((iVq|iV_q))^, and phonon GF Xph{(l,uj) = 
~{{uq\u -q))uj- The gap equation takes the form: 



-\-oo 

(k,w) = ^^ / '^^{[•^k-q - 14-q] ^tanh 



z 

2r 



+ i^(-)(c^,z|q,k-q)}[-(l/7rQ2)ImFf (q,z)]. (11) 

Here the exchange interaction Jq = 2 J (cos g^; -I- cos qy) 
induces pairing in MFA, while the Coulomb repulsion 
Vk_q suppresses the pairing. The pairing induced by 
retarded interactions is described by the kernel ([TU]). 

To estimate contributions from various interactions 
in the gap equation (jlip we consider a weak coupling 
approximation for the kernel ([TU]), i4r*^~'(aj, z|q, q') ~ 
K^~\u} = 0,z — 0|q, q'). In this approximation the 
gap equation reduces to the BCS-type form where the 
interactions are determined by the static susceptibility, 
Xq = (I/tt) X^^(dw/w)Imx(q,w): 

^C') = ^ E {'^k-q - Vk-q - |i(q)P X./(k - q) 
q 

+ |5k-q|\p/.(k - q)} |M tanh|^, (12) 



2r' 
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where = [e^ + |(p(q)|2]i/2 and ^(k) = (^2,a(k,0). 
The unimportant contribution from charge fluctuations 
Xc/(k — q) is omitted here (see later). To obtain an equa- 
tion for superconducting Tc it is sufficient to consider a 
linearized gap equation (|12p. Using a model d-wave gap 
function, ((J'(k) = A 77k, 77k = (cosfca;— cosfcj,), a linearized 
gap equation (jl2p for Tc can be written as: 



N 



-K-|<(q)rx./ + ^. 



J^SLtanh^. (13) 



The coupling constants are given by the expressions: 
" cosfc^, x7f = '^'^^'^^ 

k k 

^•^P = TfX! COSfc:!,. 



(14) 



To estimate the contribution Vc from the CI we consider 
a model for the 2D screened CI suggested in Ref. jT(i |: 



F(k) = 



1 



2^ 



(15) 



where k is the inverse screening length (|k| and k are 
measured in units of 1/a), a is the lattice constant, and Eq 
is the static dielectric constant of the lattice (in cuprates 
Eq ~ 30). For the static spin- fluctuation susceptibility 
we adopt the model as in [8|, ^] : 



X./(?,k) 



1 + [1 + (1/2) (cos k,^ + cos ky)] ' 



(16) 



Here xo = (3/4c^,)(l-<5)[(l/7V) Eq X./(q)/xo]-' is fixed 
by the condition: (Sf ) = (3/4)(l — 5) where J is a 
characteristic spin-excitation energy. The EPI coupling 
constant Vep strongly depends on the k-variation of the 
EPI matrix element |(7(k)p and a phonon dispersion in 
Xp;i(k) = l/Moj^. In particular, for a local interaction 
(/(k) = g and a dispersionless optic phonon, Wk = 
the coupling constant for the c?-wave pairing vanishes. 



ep 



0. A large electron-phonon coupling for the d- 
wave pairing can occur for a strong forward scattering, 
/c ^ in EPI (see, e.g., [i,[ll). 

Numerical integration in for the model ([T5|) gives 
for the CI coupling constant: 

K = wc 0.05 (0.11), K/^To = 0.26 (0.38), (17) 

for K = 1 (0.2) , respectively. A small ratio Vc/Vco, where 

Vco = (l/A^)Ek^(k) shows that for the d-wave pairing 
the repulsion induced by CI is remarkably suppressed. In 
particular, for Uc — 1 eV we have still a positive, though 
a small contribution from the AF exchange interaction, 
J~Vc^ 0.08 (0.02) for J = 0.13 eV. Therefore, in MFA 
we obtain only a weak coupling and a low Tc (cf. with 

Bill)- 




0.000 



FIG. 3: (Color online) Tc{S) (red solid line) compared with 
pairing induced by spin fluctuations (blue dashed line) and 
AF interaction (dotted black line) (in units of t). 



The spin-fluctuation coupling constant in (|14l) for the 
model susceptibility (IT6|) is given by, 

xTfiO = -0.66, (-0.26), xoiO = 14.8 (3.4) , (18) 

in units of 1/t = OA/us for ^ = 3.4 (^ = 1.4) at hole dop- 
ing S = 0.05 (0.30), respectively [fi]. While the spin sus- 
ceptibility Xo = Xs/(Q) Sit the AF wave vector Q is posi- 
tive and quite large, the contribution of the static suscep- 
tibility to the coupling constant XsfiO is negative 
that results in attraction mediated by spin-fluctuations 
in the equation for Tc- In the underdoped region 
with large AF correlation length ^ the spin-fluctuation 
coupling constant is quite large, while for the overdoped 
region with small ^ the coupling reduces resulting in low- 
ering of Tc- Using a conventional dispersion for elec- 
trons: t(q) = 2t{cosqx + cosqy) + At' cos qx cos qy with 
t = 0.4 eV and \t' /t\ ~ 0.2, we can estimate the spin- 
fluctuation coupling constant averaged over the Fermi 
surface, (...)f as: (1/t) {\t{q)\'^)F ~ 4i ~ 1.6 eV. Numer- 
ical estimation for the charge fluctuation susceptibility 
appears negligibly small, Xcf ^ (1/0 ^ 10"'^ which re- 
sults in a small contribution from the CI in the kernel 

m- 

The gap equation (jl2l) in strong-coupling approxima- 
tion in the imaginary Matsubara frequency w„ represen- 
tation can be written as. 



(^(k) 



-TT 



<p(q) 



[Jk-q-Vk-q-|i(q)Px./(k-q)] , (19) 



where Zk = 1 + Ak — 1 — ((i/da;)Re(I](k, tj)|i^=o is the 
quasiparticle weight. The latter is determined by the 
normal self-energy (|8]) which depends on contributions 
from all Z-channels of interactions expanded in a series 
of the Legendre polynomials Pi(cos6), contrary to the 
anomalous self-energy (|9]) where only the I = 2 channel 
contributes to the d-wave pairing. Therefore, a strong 
EPI in the / = channel resulting in a large effective mass 
renormalization, large Zk, is unimportant for the d-wave 
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FIG. 4: (Color online) 2D projection of the superconducting 
gap function <^(k) for < kx,ky < 2n. 



pairing and can only suppress the superconducting Tc 
(see also [11]). Figure [3] shows doping dependence Tc{S) 
in units of < ~ 0.4 eV for = 3 where Tc{S) induced 
by partial contributions, AF and Coulomb interactions 
in MFA oc (Jk-q — Vk-q) and spin fluctuations, cx |t(q)P 
are also shown. The maximal Tc is of the order of 150 K, 



while for Zk = 1 its value appears about five times higher. 
The gap function found for the hole concentration S — 
0.12 is shown in Fig. |4] which clearly demonstrates the 
d-wave symmetry. 

In summary, we can conclude that the superconducting 
pairing mediated by the AF exchange interaction in MFA 
is suppressed by the screened Coulomb interaction and 
only charge fluctuations cannot produce superconducting 
pairing as found in Ref. [l^]. However, spin- fluctuations 
induced by the kinematic interaction give a substantial 
contribution to the d-wave pairing and high- Tc can be 
achieved. EPI can be important for the d-wave pairing 
only for particular phonon modes having a large I = 2 
component, while polaronic effects induced by a large 
I — component of the EPI may be detrimental for su- 
perconductivity in cuprates. 
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